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NONEQUILIBRIUM RADIATION AND CHEMISTRY MODELS 
FOR AEROCAPTURE VEHICLE FLOWFIELDS 


I. Introduction 

This report covers approximately the period January 1990 thru June 1990. The primary 
tasks during this period have been the development and evaluation of various electron and 
electron-electronic energy equation models, the continued development of improved 
nonequilibrium radiation models for molecules and atoms, and the continued development 
and investigation of of precursor models and their effects. In addition! during this period 
work has been initiated to develop a vibrational model for the viscous shock layer (VSL) 
nonequilibrium chemistry blunt body engineering code. Alsoi an effort has been started 
associated with the effects of including carbon species* say from an ablator* in the 
flowfield. It should be noted that the individuals involved with these last two efforts! 
while associated with the project* are currently supported by the department. 

II. Personnel 

The staff associated with this project during the present reporting period have been 
Dr. Leland A. Carlson, Principal Investigator, and Thomas A. Gaily and Scott Stanley, 
Graduate Research Assistants. It should be noted that Mr. Gaily is currently supported by 
a NASA Graduate Student Researchers Fellowship from NASA Johnson Space Center and 
will use the results of his research on this project in his PhD dissertation. Mr. Stanley, 
who is supported entirely by this project, will use the results of his research for his 
masters thesis. In addition* two additional graduate students* Derek Green and Rajeev 
Koteshwar, have recently initiated their masters' research work in areas associated with 
the project. Currently, they are supported by departmental funds; but outside support will 
have to be found in the near future. 


III. Electron Energy Modelinc 


In the results presented in the previous progress report and in Appendix I of this 
report* the electron temperature was determined using a quasi-equilibrium free electron 
equation; and the electronic temperature was assumed to be equal to the free electron 
temperature. While it is believed that this approach is a good approximation for most 
conditions of interest in aerocapture, efforts have continued to improve the modeling of 
electron energy* and hence temperature, due to its importance in determining 
nonequilibrium ionization chemistry and radiative transfer. Specifically, studies are 
currently in progress using a combined electron-electronic energy differential equation 
which includes the effects of convection, conduction, and diffusion in addition to the 
production and loss of electron energy through elastic and inelastic collisions. The 
current full electron-electronic energy equation for the stagnation line is 
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In this equation, the viscous work terms have not been included due to the fact that they 
are of lower order. In addition, radiation effects on electron-electronic energy have been 
neglected as has electron-vibrational coupling and diffusion effects on collisions 1 energy 
exchange. The latter is expected to be small in most cases due to the rapid dissociation 
of N2 and the assumption of ambipolar diffusion. However, it might be important at some 
of the lower AFE velocities. 

It should be noted that Eqs. (E-i) is equivalent to that presented by Gnoffo (Ref. 1) 
and J. H. Lee (Ref. 2). However, it differs slightly from that presented in Ref. (3 and 4) m 
that the latter contains the additional terms 
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which arise as a result of the differences in the derivation of the species energy and 
momentum equations. It is believed that these additional terms occur as a result of using 
the more detailed approach of Chapman and Cowling (Ref. 5), and an outline of t e 
derivation is shown in Appendix II. In any event, these two terms are expected to e 
small, and their neglect in the present study should not affect the results. 


This full electron equation is included in the VSL solution sequence at the same 
location as the algebraic quasi-equilibrium electron equation was previously; and, after 
linearizing the right hand side and finite differencing along the normal coordinate, the 
resulting set of tridiagonal equations can be solved for the electron temperature, while 
this procedure is convergent, methds to improve the convergence rate are u er 
investigation. 


Results comparing temperature profiles computed with the quasi-equilibrium method 
with those obtained using the full electron equation are presented for a senes of 
velocities on Figures 1 thru 4. These results, which do not include radiation cooling or 
coupling, indicate that the primary effect of using the full electron equation is to slightly 
delay the rise in electron temperature and to move the location of peak electron 
temperature to a position slightly further from the shock front. 


One of the advantages of the new full equation formulation is that in the wall thermal 
layer the addition of the convection, conduction, and diffusion terms interact with the 
collisional terms to maintain the electron temperature equal to the translational 
temperature. In the quasi-equilibrium approach this expected behavior was not directly 
achieved but was forced by requiring Te to be less than or equal to heavy particle 
translational temperature. 


As shown on Figure 4, the full electron equation result for 8.915 km/sec is considerably 
different and cooler than the profile predicted by the quasi-equilibrium method, whi e 
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this result is still under investigation* it could be due to the low degree of ionization for 
this case and the susequent dominance of the electronic portion of the equation. In 
evaluating this result* it should be noted that the present results do not yet :nclude 
vibration electronic coupling or photochemistry effects on electron electronic temperature* 
which could change this particular profile. 

It is believed that considerable progress has been made during this reporting period in 
electron-electronic energy modeling. However* efforts will continue to improve it. In 
particular, studies are in progress to determine the importance of diffusion in collisional 
energy exchange, to assess the magnitude of photochemistry processes in the shock layer, 
to improve the solution scheme, and to include electron vibration coupling. 


IV. Noneouilibrium Molecular Radiation Models 

In the present engineering approach, nonequilibrium radiation is computed by using the 
RADICAL radiaitve analysis code and absorption coefficient model with actual species 
concentrations and correction factors to account for local thermodynamic nonequilibrium 
effects (LTNE), i.e. non-Boltzmann state populations. Previously, approximate correction 
factors for molecular radiation had been developed by Greendyke et al (Ref. 6). However, 
as alluded to in the last progress report, it is now believed that these approximate factors 
overcorrect and underestimate the actual molecular radiation. This belief is reenforced by 
the fact that experimental measurements made in molecular radiation dominated shock 
flows show a radiation intensity peak behind the shock in conjunction with the predicted 
electron temperature peak. Thus, significant depletion of the excited molecular states, as 
predicted by the theory of Ref. 6, is not expected. Consequently, it was decided to develop 
new improved molecular correction factors for molecular nonequilibrium radiation. 

After examining various approaches, a quasi-steady approach similar to that of Park 
(Ref. 7) has been developed which computes the electronic state populations associated 
with the radiating bands. Specifically, for N2, the populations of the X, A, B, a, and C 
states are computed; while for N2+ the X, A , B, and D are included. This approach has 
been developed and incoporated into the VSL code and is currently being tested. Note that 
there is no assumption concerning the existence of equilibrium between excited states and 
atoms as there was in the original factors. Thus, in this new model, both source functions 
and absorption coefficients associated with molecular band radiation will have corrections 
factors. 

The resultant molecular correction factors for the stagnation streamline of a 230 cm 
nose radius vehicle entering at 12 km/sec at 80 km are shown on Figures 5 and 6. This 
case was selected because almost the entire shock layer is in chemical nonequilibrium and 
a significant portion is in thermal nonequilibrium. Considering that unity implies no 
correction, it can be seen that for the N2 Birge-Hopfield band the correction for the 
absorption coefficient is minor but that for the corresponding source function is quite 
large in the nonequilibrium portion of the shock layer immediately behind the shock front. 
This behavior is what would "normally 11 be expected since N2(BH) involves absorption to 
the ground state. Likewise N2(i+) displays only a slight correction for the source function 
but a significant change in the absorption coefficient (from Boltzmann). This trend is also 
"expected" since N2(i+) involves the two excited states B and A. On the other hand, while 
the absorption coefficient factor for N2(2+) is similar to that for N2(i+)» the source factor 
for N2(2+) is significantly small in the chemical and thermal nonequilibrium region behind 
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the shock -front, indicating that pre-dissociation is significantly depleting the population 
of the C electronic state. 


The most interesting result, however, is that the N2+U-) radiation is only slig y 
affected by nonequilibrium phenomena, with the correction factors being essentially uni y, 
or in the case of the source function, slightly above unity. This lack of correction is in 
agreement with experiments which, at least at lower velocities, have in ica e *ivinn'i-i 
N2+U-) contribution. However, the results on Figs. 5-6 for N2+ are somewhat d 
that the number density of N2+ is only significant in the region immediately behindI the 
shock front <0.9<x/L<i.O) where the electron temperature is still increasing. Thus, most 
the N2-M1-) radiation originates from the region immediately behind the shock fron . 

Figure 7 compares, for the 12 km/sec 80 km case, the continuum radiation to thewall 
with and without the molecular nonequilibrium radiation corrections. (The dashed linei 
the result including molecular LTNE corrections.) While atomic line radiation has bee 
included in the calculation, the plot only shows the atomic and molecular continuum for 
clarity; and, as can be seen, the primary effect of nonequilibrium is to esse y 
eliminate the N2(BH) contribution between 6.5 ev and 12.77 ev, leaving primarily _o y 
atomic continuum radiation in that region. Also, nonequilibrium reduces the_N2 i-and12! ) 
bands between .75 and 4.5 ev; although, the essentially unaffected N2+d-) emission 
dominates this region. It should be noted that these trends are in agreement with the Fire 
II flight data which detected molecular radiation, particularly from the N2(l+ and - 
bands, significantly below values predicted assuming local thermodynamic equilibrium. 

Another interesting phenomena displayed on Figures 5 and 6 is that in ^e thermal 
boundary layer near the wall, several of the factors exceed unity and become large. This 
behavior indicates an overpopulation of excited states above values w ic 
predicted by a Boltzmann distribution when intuitively an equilibrium distributim would b 
expected due to the increased density near the wall. However, m actuality, the thermal 
boundary layer is in signficant nonequilibrium in that the chemical reaction rates a 
finite and cannot keep up with true local equilibrium, which le * d ® " ecies 

concentrations above equilibrium. In addition, diffusion tends to perturb the speci 
population densities and also leads to atom and ion densities above 

which in turn leads to enhanced molecular excited state popuiations. This enhance , 
however, does not lead to increased radiative emission near the wall, and, in fact, probably 
due to the low electron-electronic temperature in that region, it does not appear to affect 
the wall radiative heat transfer. 

A significant question, of course, is how important is molecular radiation at 12 Km/sec 
since originally it was believed that molecular radiation would be unimportant at t l 
condition? Radiation uncoupled stagntion point results obtained using the present metho 
with various nonequilibrium radiation assumptions are given in Table I. 

Obviously, molecular radiation at this condition, while small, comprises about fifteen 
percent of the total and should be included in predictions; and molecular LTNE effects 
slightly reduce the molecular contribution. Thus, the suggestion in the l^t report tha 
AFE conditions a conservative prediction could be obtained by including both 
atomic radiation but utilizing only atomic LTNE corrections appears to be a reasonable 

approach. 
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(At this point* it should be noted that the value listed in Table I tor the the case 
including only atomic radiation with LTNE corrections is slightly higher than values listed 
■for the same condition in the last report. Since that time* an error in the radiative 
analysis portion of the code which was introduced during the modifications to include 
LTNE effects has been discovered and corrected. It is planned to recompute the cases 
presented in the last report in the near future.) 

Results using the new electron-electronic energy and the new nonequilibrium radiation 
molecular formulations have also been obtained for the 16 Km/sec, 80 km case* and 
temperature and ionization profiles for uncoupled and radiatively coupled cases are 
presented on Figures 8 and 9. For this case, chemical and thermal nonequilibrium 
phenomena dominate the outer thirty-five percent of the shock layer; and* as can be seen 
by comparison of the two results, radiation cooling significantly cools and changes the 
chemistry in the remainder of the shock layer. Stagnation point radiative transfer 
predictions obtained with these new models are presented in Table II. 

Comparison of these results indicates that at this condition molecular radiation is 
insignificant and probably could be neglected in approximate calculations. However, 
radiation coupling and cooling is very important and needs to be included in any radiative 
heating predictions at this condition. 

Results using the new molecular correction approach have also been obtained for AFE 
CFD Points 2 and 4* and stagnation point radiative heat transfer results for various cases 
are presented in Tables III and IV. However* due to the uncertainty of the applicability of 
the full electron-electronic energy equation in its present form to these conditions* these 
results have been obtained using the quasi-equilibrium electron energy equation 
formulation. In addition* unlike the results at higher velocities* these required limitations 
on the molecular LTNE factors at the computational point adjacent to the wall. As 
mentioned elsewhere* diffusion phenomena appear to be causing excessive atom 
concentrations in that region. As a result* investigations of the present diffusion model 
in the VSL code are currently in progress. 

As can be seen, radiation cooling effects in both of these cases is small; and the 
inclusion of molecular LTNE corrections does reduce the radiative heating* which for these 
low speed conditions is primarily molecular. Finally it should be noted* that these values 
are similar to those reported in the last progress report. 


V. Noneouilibrium Atomic Radiation Models 

In the results obtained to date* local thermodynamic nonequilibrium effects (LTNE) on 
atomic radiation have been computed by applying correction factors to the absorption 
coefficient and source function values computed by the radiative analysis code RADICAL. 
The current correction factors are what should probably be termed first order 
approximations* and their form and the logic behind them is given in References 3 *6* and 8. 
Briefly, these factors assume that atomic ionization proceeds by excitation from the three 
low ground states to the high excited states followed by rapid ionization. Consequently, 
the model assumes that excitation from the ground states to the higher states is a rate 
limiting step for the ionization process and* that the excited states* because of their 
enerqy proximity to the ionized state, are in equilibrium with the free ions and electrons. 
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In contrast, Park (Ref. 7 and 9) and Kune et al (Ref. 10) handle LTNE by using a 
quasi-steady analysis in which, while rate processes between all the bound states and 
between the bound states and the ionized state are assumed finite, they are assumed to be 
fast relative to changes induced by the flowfield. Thus, at any point in a flowfield an 
equilibrium between the states will exist which is perturbed from a Boltzmann distribution 
due to radiative effects. Kune et al have performed calculations in which they specify the 
electron temperature and the total number of charged particles (defined as 2 times the 
number of atoms plus the number of ions plus the number of electrons), leaving the actual 
number of ions and free electrons to be determined as part of the unknown populations. 

Park, on the other hand, in the application of his method assumes the number of ions 
and electrons to be given by a flowfield solution (Ref. 9). Under this approach, a 
non-Boltzmann distribution can be achieved even in the absence of radiation, if the number 
of ions and electrons differs from equilibrium. To be totally correct, however, the 
excitation and ionization rates associated with each level must overall be consistent with 
the ionization rate used in the flowfield solution. 

Obviously, the present approach and those of Park and Kune et al represent the 
extremes of modeling LTNE atomic phenomena. Unfortunately, the present approach is 
overly simplified in its assumptions that the rates between the excited states and the 
free ions and electrons are infinitly fast (i.e. local equilibrium); and the detailed 
quasi—steady approaches are computationally intensive because they include a large 
number of electronic levels. In addition, the latter are sensitive to the choice of the 
individual rates; and it is difficult to know which rate to ad^jst when comparing with 
experimental results and attempting to improve the correlation. Consequently, work has 
been initiated during this reporting period on an improved LTNE atomic model. 

After extensively reviewing the work on argon of Foley and Clarke (Ref. 11), Nelson 
(Ref. 12), etc. and the air and nitrogen work of Park (Ref. 7), Kune and Soon (Ref. 10), and 
others, it has been decided to model high temperature nitrogen by subdividing atomic 
nitrogen into two species. The first, termed Ng, for N ground, represents the nitrogen 
atoms in the first three low lying electronic states of nitrogen. The second, termed N* or 
N excited, representes those nitrogen atoms populating the remaining upper electronic 
states. The relative densities of these subspecies will then be determined by appropriate 
reaction rates between themselves, N+, e-, etc. It is believed that this approach will be a 
significant improvment over the present model in that it will allow a finite rate of 
ionization from excited states while retaining the fundamental two step ionization 
process. In addition, by determining the excited state number densities directly from the 
flowfield computation, the appropriate atomic LTNE factors should be "directly 
obtainable and more accurate. 

Initially, it is planned to use the collisional reaction rate system shown in Table V. In 
general, reaction rates for the first seven reactions are well known. However, the rates 
for the electron-atom excitation and electron-atom ionization reactions, numbers 8-10, 
need to be determined. Currently, atom-atom excitation and photo-excitation 
photo-ionization are not included since it is believed that these reactions are of second 
order in the stagnation region. However, it is planned to include them later. 

In this system, care must be taken to properly formulate the species enthalpy of Ng and 
N*. Specifically, 
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where 
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For equilibrium conditions* these expressions reduce to the proper forms where 


h - £ & k 
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As mentioned above, effective reaction rates have to be obtained for reactions (8) - 
(10). While in principle, these could be extracted from the work of Park (Ref. 13), the work 
in Ref. 10) appears to contain information based upon more recent data. Furthermore, it 
appears to yield excitation rates more compatible with relaxation data behind shock waves. 
Consequently, a method has been developed and a computer program written to determine 
from the detailed data of Ref. 10, effective forward rates for reactions (8)-(10). While 
complete details of the method and results will be presented later, a preliminary set of 
results is presented in Figure 10. 

Also shown on Figure 10 is the rate of Wilson successfully used in Ref. 3 in conjunction 
with the first order LTNE model. As can be seen the preliminary rates for the new model 
faster for excitation from the ground state but are finite for ionization from the 
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excited state to the continuum. Thus» they appear to have the right trend and magnitude. 
In addition* included on the figure is the effective ionization rate from the ground state 
directly to the continuum. As previously postulated, this rate is considerably slower than 
the excitation rate. Finally# the ground to excited forward rate is about two orders of 
magnitude slower than that which it is believed would be obtained from using the detailed 
rates in Ref. 13. 

Currently, various methods of curve fitting these rates and incorporating them into the 
radiation coupled nonequilibrium VSL code are being investigated. It is anticipated that 
the results of these studies will be reported in the next progress report. 

Once the chemistry model involving excited species has been developed, the next step 
is to determine the appropriate LTNE factors which should be utilized in the radiative 
analysis code. Considerable progress has been made in obtaining these factors during this 
period. However, since they are still under development and being “checked 11 , the results 
which follow should be considered preliminary. The derivation of these factors will be 
presented at a later date. 

For continuum processes involving absorption by an excited state, the absorption 
coefficient factor is 



and the factor on the source function is . 
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Similarly, for continuum absorption involving a "ground" state, the absorption LTNE factor 
is 



and the corresponding source function factor is 
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For line processes involving absorption into an excited state* the present theory yields an 
absorption LTNE factor of 
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while -for this case the source function is unchanged. On the other hand* if the line process 
involves absorption into one of the ground states* the absorption factor is 
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and the source function also has^factor 
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It should be noted that if the N* species is in equilibrium with N+ and e- and the number 
density of Ng is assumed to be that of N, these expressions reduce to those used with the 
first order model. 


While the effects of using this new chemistry model and improved LTNE radiation model 
will not be Known until they have been incorporated into the VSL code and extensively 
tested and compared to previous resultsi some preliminary studies of the effect of the 
LTNE factors have been conducted using temperatures and species densities previously 
obtained. These results indicate that the absorption factors are often larger than those 
obtained with the first order model. Since the new model has a faster excitation rate 
between ground and excited states and a finite rate between excited levels and the 
continuum* it should lead to higher population densities in the excited states. Thus* the 
preliminary results appear to have the correct trend. Consequently* it is quite possible 
that radiative heat transfer values obtained with the new LTNE model may be higher than 
those previously predicted. However* since these preliminary studies do not include the 
coupling effects on chemistry and flowfield profiles which different values of LTNE 
factors will induce* estimates of the magnitude of the change would be purely speculative 
at this time. 


VI. Precursor Studies 


As noted in the last progress report* precursor phenomena leads to the appearance 
ahead of the vehicle shock front of ionized species, free electrons, and thermal 
nonequilibrium each of which possibly might have an affect on the subsequent post-shock 
nonequilibrium phenomena. During this reporting period the effort to develop a good 
engineering model for precursor phenomena in nitrogen has continued and has led to 
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changes in the absorption coefficient model, improvements in the soiution scheme, the 
development of appropriate radiation view factors, and the inclusion of elect 
nonequilibrium. 

Radiation View Factors 

In computing shock layer radiative phenomena, the usual approach in peering 
calculations is to use the tangent slab approximation. Since the ratio of sho-klay 
thickness to vehicle radius or diameter is small, this approach is appropriate. * * 

the precursor region ahead of the shock front, important phenomena occur at significan 
distances from the vehicle; and at those points the radiating shock layer ^ the vehicle 
only comprises a small portion of the spherical field of view. In other words, asthe point 
of interest in the precursor moves away from the shock front, the shock and y 
appear to be infinite slabs; and the actual solid angle over which the radiation should b 
spacially integrated must be properly computed. 

By assuming that there is no emission in the precursor, it can be shown that the 
approriate radiation expression for a point in the precursor is 

\ r ^ (Tv seed) 1 


-217 (fs&irrt*)* 


£ 3 ir y -u) -tiiv J 


(P-1) 


where & is one-half of the angle subtended by the body as viewed from the pQint in the 
precursor. The complete derivation of this expression will be given in detail at a later 

date. 

It should be noted that Eq. (P-1) is essentially the tangent slab expression except that 
each term has been modified or corrected by a view factor depending upon vehicle size, 
location of the point of interest, and frequency. The typical magnitude and variati 
these factors as a function of distance in front of the shock is shown on Figures ii and 12. 
On this figure, the curve denoted “wall correction” refers to the view factor on the first 
term of Eq. (P-i) while that marked "shock layer correction” is the factor on ^e second 
term. Since it is anticipated that the radiation from the shock layer, rather\ h * ^ 
wall, dominates the precursor phenomena, the present engineering precursor mo 
the shock layer view factor on all terms as an approximation, i.e. 
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where AF^ is the view factor and cu is the total radiative flux at a point assuming a 
tangent slab. ” 

In the energy equation, the term involving radiation appears as a divergence o-f the tlux 
and is defined to account for the net absorption-emission at a point. However» simple 
differentiation of Eq. (P-2) would yield 
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In this expression» the first term is the change in the radiative flux due to absorption or 
emission, but the second term is the change due to geometry and should not be included in 
the energy equation. If it were included, an essentially transparent radiation would 
appear to be absorbed due to the spatial variation in the view factor. (Consider the case 
where is constant.) 

Absorption Coefficients 

Durinq this reporting period, there have only been a few minor changes in the precursor 
absorption coefficient model. Originally, predissociation from both the 1st and 2nd 
positive bands of N2 was permitted since it was theoretically possible. However, 
subsequent review of the literature indicates that in the present temperature range of 
interest in the precursor that predissociation is involved in the the photodissociation of 
nitrogen only from the state following absorption (photoexcitation) in the 

Lyman-Birge Hopfield Band. Therefore, predissociation is no longer included through the 
N2 1st and 2nd positive bands. 

The second change is that the Vegard-Kaplan band is no longer included in the radiative 
calculations because, first, there is little evidence in the literature of its existence and, 
second, numerical studies conducted as part of this research have shown that is has 
negligible effect on precursor phenomena. 

It is recognized that while these processes appear to be unimportant, they could have 
been retained in the model for completeness. However, their inclusion does require a 
measureable amount of computer time in any given case; and thus, in the interest of 
efficiency, they have been dropped. 

Consequently, the following radiative processes are currently included in the precursor 
model: 
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Process 


Frequency Range 
(ev) 


N Low Frequency Ionization 
N High Frequency Ionization 
Free - Free 
N2 Ionization 
N2+ <1-) 

N2<B i rge-Hopf i e 1 d) 

N2(1 + > 

N2(2+) 

N2 (LBH) 


h v 
h v 
h v 
h v 

2.23 < h v 
6.5 < h v 
.75 < h v 
.75 < h v 
h u 


> 0 

> 10.8 ev 

> 0 ev 

> 8.24 eu 

< 4.46 eu 

< 12.77 ev 

< 4.5 ev 

< 4.5 ev 

> 4.77 ev 


Solution Scheme 


Originallyi the solution scheme -for the precursor region solved the governing 
differential equations assuming a constant spatial step size. However, this approach in 
many cases either failed to provide sufficient detail in regions of high gradients or led to 
excessive memory and CPU requirements. Consequently, the solution scheme has been 
modified to allow a variable step size that is determined by changes in the flowfield 
properties. Typically, the step size is limited to a maximum value corresponding to a 
maximum of 15% percent change in any flowfield property; while the minimum value is 
selected to yield at least a 5% change. This approach significantly reduces solution times 
while placing large numbers of points in regions of high gradients. 


Electron Thermal Noneauilibrium 

The primary objective of the precursor research during this reporting period has been 
to develop and include in the nitrogen precursor model an appropriate model for electron 
thermal nonequilibrium. Such a model is important because the precursor radiative 
phenomena are strongly determined by the electron temperature. 

While in the shock layer it is frequently possible to use a free electron temperature 
and assume that due to collisional phenomena that the electronic temperature is equivalent 
to the electron temperature, such an approach in the precursor is tenuous due to the 
overall low density and low number of free electrons. While theoretically a three 
temperature model, T, Te, and Telec, could be conceptualized, the electron electronic 
energy exchange expressions are not well known or understood. Consequently, it has been 
decided to utilize in the precursor model, two temperatures, T and Te=Telec, and a 
combined electron-electronic energy equation. 

In the development of a combined electron-electronic energy equation, particularly 
considering that the dominant processes in the precursor are radiative, the manner in 
which each radiative process affects electron energy must be considered. Specifically, 
free-free processes affect the free electron translational energy, while in bound-free 
photoionization the energy involved is divided between that required for ionization and 
that associated with the translational energy of the created electron. Similarly, atomic 
line phenomena only affects the electronic energy, but the energy of molecular band 
absorption is distributed between electronic, vibrational, and rotational energy changes. 
Finally, bound-free photodissociation absorption involves the energy of dissociation plus 
the kinetic energy of the created atoms. Obviously, the inclusion of the appropriate 
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portions o i each of these processes in the electron-electronic energy equation is 
somewhat subtle. 

During this period* two approaches to the electron electronic energy equation have been 
investigated* and in both cases thermal conduction* electron vibrational coupling* and 
diffusion effects have been neglected. The first approach is the “usual 11 technique of 
combining the free electron and electronic energy equations. By defining a new variable as 



(P-4) 


the combined equation can be written as 

* (pe“') + &.(pe:“i)= - fif.■« + X ^ fQ <- 

^ -v £ <p-5) 

1 - 



Here the last term in the first equation accounts for the radiative energy that goes into 
ionization in photoionization processes, while the second term in the last equation 
accounts for electronic energy changes due to molecular band transitions. In principle, the 
free-free* bound-free photoionization, and atomic line flux divergences could be 
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individually obtained from the radiation model* RADICAL. However* the last integral still 
requires the computation of total radiative flux. 


The second approach is to directly utilize the total radiative flux computed by the 
modified version of RADICAL and eliminate from it those portions which do not directly 
affect electron or electronic energy. An appropriate energy equation for this case can be 
obtained by multiplying each species continuity equation by the species zero point energy* 
sumino over all species and adding the resulting equation to the combined free electron 
electronic equation. Then by defining* 
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a combined electron-electronic energy equation is ( >„ 
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Notice that the last four terms essentially correct the divergence of the radiative flux for 
those radiative portions which do not affect the electron or electronic energy. By 
properly examining each of the radiative flux terms* Equation (P-7) can be rewritten into 


the form 
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Details of the derivation of these equations will be presented in a future report. 


With respect to this electron-electronic energy equation and model the following points 
should be noted. First* the second approach decreases the number of required frequency 
inteqrations and modifications to RADICAL. Second* these equations are also applicable to 
the shock layer if it is desired to include in the shock layer the effects of photoprocesses 
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on electron electronic energy. Finally* it is believed that the above equation and model* 
when combined with our modified version of RADICAL* represents the most complete 
electron electronic energy and precursor model to date for a complicated diatomic gas such 
as nitrogen. 

Test Case 


In order to develop the precursor model* it has been decided to use as a test condition 
the shock layer at 12 km/sec for a 230 cm nose radius body at 80 km altitude in nitrogen. 
In the last progress report* results were presented which neglected radiative cooling in 
the shock layer* assumed local thermodynamic equilibrium in the shock layer* and assumed 
thermal equilibrium in the precursor (i.e. Te = T). During this reporting period* preliminary 
results have been obtained for this test condition. However, these new results include 
thermal nonequilibrium* radiative cooling* atomic line* atomic continuum* and molecular 
radiation in the shock layer. They also include local thermodynamic nonequilibrium 
phenomena for the atomic radiative processes in the shock layer. In the precursor, these 
results include radiation view factors* thermal nonequilibrium (Te not equal to T)* 
photodissociation and photoionization chemistry, and continuum emission and absorption 
processes. However* the precursor results do not yet include line effects or collisional 
chemistry. 

Initially, the populations of excited states in the precursor were assumed to be 
determined by a Boltzman distribution at the electron temperature. While not strictly 
correct at the low densities being considered, it was believed that such an approach would 
be a reasonable approximation and not lead to any significant errors* particularly 
considering the low electron temperatures expected in the precursor zone. However, 
examination of the results indicated that the electron temperatures were sufficiently high 
to lead to emission in the region immediately in front of the shock, with a subsequent 
decrease in electron temperature. 

Consequently, a collision limiting correction (Ref. 14) has been included in the program 
to more realistically predict the populations of the excited states of the molecules; and 
some preliminary results are shown on Figures 13 thru 18. Figure 13 shows the heavy 
particle temperature variation in the precursor zone; and in contrast to the results shown 
in the last progress report assuming thermal equilibrium, the present results indicate that 
heavy particle temperature is essentially constant in the precursor zone. This result is 
expected and is in agreement with multi-temperature calculations for argon (Ref. 11). 

Likewise* the variations in pressure and density* portrayed on Figures 14 and 15* also 
are* for the test case, essentially constant in the precursor region. In this case, the 
slight increase in pressure ahead of the shock front is less than shown in the last 
progress report because the previous results did not include cooling or local 
thermodynamic nonequilibrium effects in the shock layer. Thus* the previous precursor 
results were for a case with a more strongly emitting shock layer* which induced more 
photoionization and dissociation. In the present case* the inclusion of cooling and 
radiative nonequilibrium effects in the shock layer results significantly reduces the 
radiation to the precursor zone. 

The variation of the electron-electronic energy in the precursor is shown on Figure 16. 
As can be seen, as the flow approaches the shock front, it absorbs radiation from the bow 
shock layer; and the electron-electronic energy increases exponentially. Similarly, the 
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species concentrations, shown on Figure 17, also vary exponentially with distance from the 
shock front. Notice that while low, there is measurable dissociation at distances greater 
than 150 shock layer thicknesses, indicating that photodissocition of N2 occurs primarily in 
the tar precursor. In addition, there is significant ionization, on the order of IE-04 and 
primarily N2+, in the near precursor immediately in front of the shock wave. Althoug e 
effect of these pre-shock ions and electrons on the post-shock flowfield has yet to be 
determined, it should be noted that the present results are for » irelatively cooJ| shoe 
layer (10,000 K) and yet they indicate that immediately behind the shock front, an e 
density on the order of 4E11 particles/cubic cm could be expected. Such a level might 
affect both chemical and radiative nonequilibrium phenomena in the nonequilibrium cone 
behind the shock front. 

Finally, the electron temperature variation in the precursor is shown on Figure IS. In 
examining this result, it should be remembered that electron temperature is a measureof 
the average kinetic energy associated with each electron at a given location. Thus, the 
high temperature predicted for the far precursor indicates that electrons seated via 
photoionization far away from the vehicle were created by high energy radiation. Howeve , 
as shown on Figure 19, the number of electrons in the far precursor is extremely small. 
Then, as the flow approaches the vehicle the electron temperature slowly increa ** s t0 4 
value above 4000 K until about 50 shock thickness, where it begins to decrease slightly. 


Originally, it was suspected that this electron temperature decrease was due to 
emission from the near precursor created by the assumption that the electronic states o 
N2+ etc. were populated according to a Boltzmann distribution. This assumption predicts 
higher populations for the excited states than can be realistically maintained by 
collisional excitation at the freestream densities and leads to ehanced emission. This 
supposition has been partially verified by the introduction of the collision limiting mo e , 
which predictes electron temperatures (shown on Fig. 18) near the shock front about 
higher than those obtained without it. Interestingly, electron temperature is the on y 
quantity measureably affected by the introduction of the collision limiting model. 


Currently, it is believed that the electron temperature decrease shown on Figure 18 is 
due tD the exponential growth of the number of electrons in the near precursor region 
combined with the fact that these electrons are created by photoionization involving 
radiation near and slightly above the ionization threshold. Consequently, the coated 
electrons have "low" kinetic energy compared to those formed in the far precursor; and the 
average free electron kinetic energy or electron temperature decreases. However, when 
atomic line radiation is included in the precursor model, it is anticipated that the electron 
temperature will increase somewhat in the near precursor due to the absorption of energy 
into the electron bands. On the other hand, the inclusion of collision chemistry in the near 
precursor would lead to some ionization by electron impact, which in turn would cause a 
decrease in electron temperature. Thus, the present profiles may be reasonably realistic 
for a nitrogen freestream. In any event, these results do indicate that enhanced elec ron 
temperatures on the order of 4000 K and measurable ionization levels due to N2+ do exist 
in *the reoion immediately in front of the shock layer. 


Future Plans 

During the next reporting period it is planned to obtain further precursor solutions, 
probably for the velocity range of 12 to 16 km/sec and for altitudes of 70 to 80 km. It is 
anticipated that precursor phenomena will be significantly increased as altitudes are 
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decreased and velocities are increased. Based upon these results# a series of parametric 
studies will be conducted on the stagnation region shock layer using the nonequilibrium 
radiation coupled multi-temperature VSL code in order to determine the effect of these 
precursors on the shock layer chemistry# flowfield profiles# and# in particular# the 
radiative heat transfer to the body. In addition# since the present precursor 
electron-electronic energy model has been formulated to include line radiation# 
consideration will be given to including atomic lines and determining their effect on the 
precursor. However# it is not anticipated that collisional chemistry will be included in the 
precursor model during the next reporting period. 

VII. Publications 


While no new publications associated with the project were issued during this reporting 
period# AIAA Paper 89-1729, "The Effect of Electron Temperature and Impact Ionization on 
Martian Return AOTV Flowfields," by L. A. Carlson and T. A. Gaily # has been accepted for 
publication in the Journal of Thermophysics and Heat Transfer . It is anticipated that this 
article will appear during the latter part of 1990. 

In addition# an abstract of a proposed paper entitled "Nonequilibrium Chemical and 
Radiation Coupling Phenomena in AOTV Flowfields" has been submitted and accepted for 
presentation at the 29th Aerospace Sciences Meeting to be held in Reno# Nevada in January 
1991. A copy of this abstract is included in this report as Appendix I for your information 
and approval. 

During the next reporting period# it is anticipated that one or more abstracts will be 
submitted for possible paper presentations at the AIAA Fluid and Plasmadynamics and 
Thermophyics Conferences to be held in June 1991 in Honolulu# Hawaii. 


Table 1 — Stagnation Point Radiative Heat Transfer for Test Case I 
(U = 12 km/sec, Altitude = 80 km, Nose Radius = 2.3 m) 


Atomic Radiation and Atomic LTNE Corrections Only - 14.0 watts/sq cm. 

Atomic and Molecular Radiation 

with Atomic LTNE Corrections Only 18.1 watts/sq cm. 

Atomic and Molecular Radiation with 

Both Atomic and Molecular LTNE Corrections — 1<4.5 watts/sq cm 
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Table II — Stagnation Point Radiative Heat Transfer -for Test Case III 
U = 16 km/sec, Altitude = 80 km, Nose Radius = 2.3 m 


Atomic Radiation and Atomic LTNE Corrections Only — 

No Radiation Gasdynamic Coupling 

331 watts/sq cm. 

Atomic and Molecular Radiation 
with LTNE Corrections or. Both 

No Radiation Gasdynamic Coupling 

— 

340 watts/sq cm. 

Atomic and Molecular Radiation with 

Both Atomic and Molecular LTNE Corrections 

Radiation Gasdynamic Coupling Included 

189 watts/sq cm 


Table III— Stagnation Point Radiative 

U = 8.915 km/sec, Altitude = 77.9 km, 

Heat Transfer -for AFE CFD Point 2 
Nose Radius = 2.3 m 

Atomic and Molecular Radiation 
with Atomic LTNE Corrections Only 

No Radiation Gasdynamic Coupling 

— 7.0 

watts/sq cm. 

Atomic and Molecular Radiation 
with LTNE Corrections on Both 

No Radiation Gasdynamic Coupling 

— 1.7 

watts/sq cm. 

Atomic and Molecular Radiation 
with Atomic LTNE Corrections Only 
Radiation Gasdynamic Coupling Included 

— 5.7 

watts/sq cm. 

Atomic and Molecular Radiation with 

Both Atomic and Molecular LTNE Corrections " 1*7 

Radiation Gasdynamic Coupling Included 

watts/sq cm 
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— 

Table IV — Stagnation Point Radiative Heat 
<U = 9.326 km/sec, Altitude = 75.15 km, 

Transfer -for AFE CFD 
Nose Radius = 2.3 m) 

Point A 

— 

Atomic and Molecular Radiation 
with Atomic LTNE Corrections Only 

No Radiation Gasdynamic Coupling 

— 8.8 watts/sq 

cm. 


Atomic and Molecular Radiation 
with LTNE Corrections on Both 

No Radiation Gasdynamic Coupling 

3.2 watts/sq 

cm. 


Atomic and Molecular Radiation 
with Atomic LTNE Corrections Only 

Radiation Gasdynamic Coupling Included 

7.5 watts/sq 

cm. 

_ 

Atomic and Molecular Radiation with 

Both Atomic and Molecular LTNE Corrections 
Radiation Gasdynamic Coupling Included 

— 3.2 watts/sq 

cm 
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Fig, 1 — Stagnation Line Heavy Particle and Electron Temperature Profiles Using the 
Complete Electron-Electronic Energy Equation and the Quasi-Equilibrium Electron Energy 

Equation 

U = 16 km/sec, H = 80 km, Rnose = 2.3 m 
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Fig. 2 — Stagnation Line Heavy Particle and Electron Temperature Pro-files Using the 
Complete Electron-Electronic Energy Equation and the Quasi-Equilibrium Electron Energy 

Equation 

U = 14 km/sec, H = 80 km, Rnose = 2.3 m 
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Fig. 3 — Stagnation Line Heavy Particle and Electron Temperature Protiles Using the 
Complete Electron-Electronic Energy Equation and the Quasi-Equilibrium Electron Energy 

Equation 

U = 12 km/sec, H = 80 km, Rnose = 2.3 m 
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Fig. 4 — Stagnation Line Heavy Particle and Electron Temperature Pro-files Using the 
Complete Electron-Electronic Energy Equation and the Quasi-Equilibrium Electron Energy 

Equation 

U = 8.915 km/sec, H = 77.9 km» Rnose = 2.3 m 


































and Te (deg K) *10* 



Fig. 8 - Stagnation Line Heavy Particle and Electron-Electronic Temperatures 
With and Without Radiation Gasdynamic Coupling 
U = 16 km/sec, H = SO Km, Rnose = 2.3 m 
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Temperature Variation in Precursor 
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- Heavy Particle Temperature Variation in the Precursor 
U = 12 km/sec, H = 80 km, Rnose = 2.3 m 
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Fig. 14 - Pressure Variation in the Precursor 
U = 12 km/sec, H = 80 km, Rnose = 2.3 m 









Density Variation in Precursor 


V = 12 Km/ 
Alt = 80 Km 


51 101 

X/Xshock ( —) 


151 


Fig. 15 — Density Variation in the Precursor 
U = 12 Km/sect H = 80 Km, Rnose = 2.3 m 















Species Concentrations in Precursor 



Fig. 17 — Species Concentration Pro-files -for the Precursor 
U * 12 km/sec, H = 80 km» Rnose = 2.3 m 
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Fig. 18 — Electron-Electronic Temperature Variation in the Precursor 
U = 12 km/sec, H = 80 km, Rnose = 2.3 m 
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NONEQUILIBRIUM CHEMICAL AND RADIATION COUPLING 
PHENOMENA IN AOTV FLOWFIELDS 

AN EXTENDED ABSTRACT 


Leland A. Carlson* and Thomas A. Gaily** 
Texas A&M University 
College Station, Texas 


SUMMARY 

In order to investigate radiative and chemical nonequilibrium effects, flowfield solu¬ 
tions for a wide range of AOTV flight conditions have been obtained. These solutions 
have been calculated using a viscous shock layer method which includes the effects of 
chemical and radiative nonequilibrium, thermal nonequilibrium, viscosity, heat conduc¬ 
tion, diffusion, and radiative-gasdynamic coupling. The effects on radiative heat transfer 
and flowfield properties due to radiative coupling and local thermodynamic nonequilibrium 
effects are shown and discussed. The variation of radiative heating rates with velocity and 
altitude are also presented for the different coupling and radiative nonequilibrium models 
investigated. 


INTRODUCTION 

In the future, various space programs will be conducted which will require the efficient 
return of large payloads to low earth orbit from missions to the moon or to planets like 
Mars. To accomplish this task, the return vehicles will utilize aerocapture techniques 
that will involve reentry and deceleration at high altitudes; and in order to design these 
vehicles, a thorough understanding of the physical phenomena will be required. Because of 
the high altitudes associated with aerocapture, the vehicle flowfields will be dominated by 
chemical, thermal, and radiative nonequilibrium phenomena. Thus, the primary purpose 
of the present study is to develop an engineering flowfield model suitable for high altitude 
AOTV flowfields having extensive chemical, thermal, and radiative nonequilibrium and to 
use this model to investigate the magnitude, extent and coupling between these phenomena. 


* Professor, Aerospace Engineering Dept. 

** NASA Graduate Research Fellow, Aerospace Engineering Dept. 
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METHODS 


The flowfield model used in this investigation is a viscous shock layer analysis which 
includes the effects of chemical nonequilibrium, multi-temperature nonequilibrium (elec¬ 
tron and heavy particle), viscosity, heat conduction, diffusion, and radiative-gasdynamic 
coupling. The basic code is the same as used for the previous study presented in Ref. 
(1), but a number of additions and modifications have been made since the earlier study. 
First, the VSL code has been coupled with the radiation analysis routines of the NASA 
Langley program, RADICAL^, which is described below. Second, the chemical reaction 
rate input data has been changed to allow the use of a single reaction rate, kf or k^, and 
the equilibrium constant, Keq, rather than using both the forward and backward rates. 
This modification was deemed necessary since the ratio of the experimentally determined 
rates k£ and k^ often deviates significantly from the theoretical Keq- With this modifica¬ 
tion the species concentrations in the near equilibrium regions of the flowfield are now in 
agreement with the results predicted from strictly equilibrium analysis. Third, the effects 
of multi-temperatures on the shock jump conditions and thermodynamic state variables 
have been improved from those at the time of Ref. (1). 

As mentioned, the radiation analysis package from the program RADICAL has been 
coupled to the viscous shock layer flowfield, giving the ability to calculate flowfield solu¬ 
tions with the effects of radiative cooling present. The radiation analysis in RADICAL is 
a detailed method which includes atomic continuum radiation, molecular band radiation, 
and atomic line radiation for the standard CHON (carbon, hydrogen, oxygen, nitrogen) gas 
system. The original method uses the species number densities and assumes a Boltzmann 
state distribution to calculate the excited state number densities for each species, and, 
from this data, the individual radiative absorption coefficients. Using the theory of excited 
state depletion under nonequilibrium conditions presented in Ref. (3) and Ref. (4), the 
present radiation analysis uses atomic local thermodynamic nonequilibrium (LTNE) radi¬ 
ation correction factors which account for state population deviations from that predicted 
by a Boltzmann distribution. A similar set of molecular corrections factors can be deduced 
from this theory, but the authors question the validity of applying the theory to molecular 
vibrational states, and an alternate method of determining molecular state populations in 
under development (see EXPECTED RESULTS). 

More details of the flowfield and radiative methods will be presented in the final paper, 
including equations, relevant figures, and tables. 




TYPICAL RESULTS 


In this section results obtained using the nonequilibrium coupled radiation viscous 
shock layer model are presented. At this stage, these results are being primarily used to 
gain insight into the phenomena affecting the flowfield chemistry and radiation. In all 
cases, results axe for the stagnation line, utilize fifty-one points between the shock and 
the wall, and the freestream is N 2 . The nonequilibrium chemistry model is the Case II 
Nitrogen Reaction System presented in Ref. 1, with the reverse rates being obtained via 
the equilibrium constants. In addition, the wall has been assumed to be a noncatalytic 
black body at 1650 deg K. This latter assumption has been used for convenience even 
though it is recognized that for many of the cases of interest that the heat transfer load 
will be more than adequate to induce ablation and to raise the wadi temperature to values 
several thousand degrees higher. 

AFE CFDPoint 2 

This condition corresponds to what is often referred to as the “max Q” computational 
point for one of the initial AFE trajectories at which the freestream velocity is 8.915 
km/sec, freestream pressure is 15.715 dyne/cm^ and temperature is 197.101 K. For this 
case, the nose radius has been assumed to be 2.3 meters, and both atoms and molecules 
have been included in the radiation calculations. Heating results are presented in Table I 
and the stagnation temperature and composition profiles are shown as Figure 1. 

As can be seen, the results include electron thermal nonequilibrium but they assume 
vibrational equilibrium. Since for the AFE vibrational nonequilibrium effects will be sig¬ 
nificant and may affect the chemistry and the radiation, it is planned to include in the 
near future a vibrational nonequilibrium model in the nonequilibrium radiating VSL code. 
However, the present results assume Tyj^ = 

As shown on Figure 1, the electron temperature rapidly rises behind the shock front 
and equilibrates with the heavy particle temperature. However, as evidenced by the con¬ 
tinual decrease in temperature and the variations in composition across the shock layer, 
the stagnation flow for this case is always in chemical nonequilibrium. Also, the wall ther¬ 
mal layer comprises approximately twenty percent of the 10.8 cm thick shock layer. It 
should be noted that the electron temperature and composition results shown on Figure 1 
axe very similar to unpublished results previously obtained for this case using the inviscid 
AFETE code (Ref. 5) and an earlier version of the VSL code (Ref. 1). 

As would be expected for this case, the radiation coupling effects for this case are 
very small and cannot be detected on plots comparing uncoupled and coupled solutions. 
However, as shown on Table I, there is a small amount of cooling, which can be observed 







by comparing corresponding cases such as uncoupled corrected with coupled corrected etc. 
Also, as mentioned above, the radiative heating results for this case include molecules 
without any LTNE factors; and thus they should be conservative with respect to the 
radiative heating estimates, which range from 4.76 to 7.21 watts/cm 2 . Note that the 
usage of the atomic LTNE corrections factors only reduces the radiative heating about 0.7 
watts/cm 2 , indicating that most of the radiation is probably molecular. 

Since results have been obtained for this case previously using different versions of 
the method, the sensitivity of the results to various factors can be observed. For example, 
with the original VSL code, forward and reverse rate reaction chemistry, and the eight step 
radiation model used in Ref. (1), the non-radiation coupled results were with no LTNE 
corrections 29.7 watts/cm 2 and 22.4 watts/cm 2 with atomic LTNE corrections only. For 
these cases the electron temperature profile in the nonequilibrium zone was slightly higher 
than those shown on Figure 1 due to the influence of different reaction rates. On the 
other hand, after shock jump and enthalpy changes involving electron temperature were 
incorporated but before the forward and reverse rates were replaced with the forward 
rate and an equilibrium constant, the results for columns (3) and (4) of Table I were, 
using RADICAL, 15.9 and 12.6 watts/cm 2 respectively. Thus, radiative heating results 
are sensitive to small details in the methods, the radiation model, and in particular to 
the chemistry and electron temperature. It should be noted that the change from forward 
and reverse rates to only a forward rate and a computed equilibrium constant significantly 
changed the equilibrium temperature and composition. 

Since the results shown on Figure 1 and Table 1 for this case are in better agreement 
with the detailed inviscid results obtained using AFETE than previous VSL values and 
since they have been obtained using better models, they axe probably reasonable estimates. 
However, the effects of vibrational nonequilibrium and chemical nonequilibrium on molec¬ 
ular radiation have yet to be included. Nevertheless, since the latter effects should most 
likely reduce radiative heating, the presents estimates may be conservative. 

AFE CFD Point 4 

It is believed that this condition corresponds to a “max Q” computational point for an 
AFE trajectory associated with a heavier vehicle at which the freestream velocity is 9.326 
km/sec, freestream pressure is 26.4 dyne/cm 2 and temperature is 200 K. For this case, 
the nose radius has also been assumed to be 2.3 meters, and both atoms and molecules 
have been included in the radiation calculations. Heating results are presented in Table 
I and the stagnation temperature and composition profiles are shown Figure 2, for which 
the shock layer thickness is 10.5 cm. 





As can be seen, these results are very similar to those for the CFD Point 2 case, and 
the remarks concerning that case apply equally well here. In general, the predicted heating 
rates, which do not include molecular nonequilibrium radiation factors, are higher than 
those for CFD Two. 

U = 14.5 km/sec, Altitude = 65 km 

In a recent paper (Ref. 6) results have been presented for the stagnation line of a 
one meter nose radius body at a trajectory point possibly representative of an earth entry 
return from Mars. These results include chemical nonequilibrium, thermal nonequilibrium 
assuming that vibrational, electronic, and electron nonequilibrium can be represented by a 
single temperature, and uncoupled nonequilibrium radiation computed by a detailed model 
that includes the molecular continuum and atomic lines. 

Using this model, the investigators obtained for this trajectory point an uncoupled 
radiative heating rate of 1700 watts/cm 2 , a shock standoff distance of 5.7 cm, and a 
post-shock chemical nonequilibrium zone 1.1 cm thick. In this nonequilibrium zone, the 
electron-vibrational temperature never significantly exceed the equilibrium temperature. 
They also stated that most of the radiative heating was from the ultra-violet below 2000 
A, that it originated from the nonequilibrium region behind the shock wave, and that very 
little was absorbed in the wall thermal layer. The latter is different from previous beliefs by 
some researchers (Ref. 7), but it is in agreement with the approximate studies of Carlson 

(Ref. 3). 

In addition, in Ref. 6 comparisons were made with results obtained using the RASLE 
code (Ref. 8), which is an equilibrium viscous shock layer code using a radiation model also 
based upon RADICAL. Using the RASLE code, Ref. 6 obtained for the same case a shock 
standoff distance of 3.5 cm and a radiative heating rate of 970 watts/cm 2 . The authors 
attributed the differences to nonequilibrium chemistry effects and the RASLE radiation 
model, asserting that the latter smeared atomic lines and therefore obtained incorrect 
results. 

Considering these discrepancies, it is believed that it would be valuable to apply the 
present model to this trajectory point; and temperature and composition profile results 
for the case including radiative cooling are presented in Figure 3. For this case, the shock 
standoff distance is 3.4 cm; and, as can be seen, most of the shock layer is in chemical 
equilibrium. 

The difference in shock standoff distance between the present results and those of Ref. 
6 is believed to be primarily due to the electron temperature profile and its subsequent 
effect on chemistry. In Ref. 6, T e is low in the region behind the shock front, possibly due 




to combining electron-electronic with vibrational phenomena. However, as shown in the 
composition profiles, diatomic species are insignificant over most of the shock layer and 
ionization dominates the chemistry. Thus, in the present case, an electron temperature 
model which is strongly influenced by collisional and ionization phenomena is used. Figure 
3 shows that the result of using such an approach is an electron temperature which in the 
nonequilibrium zone behind the shock front significantly exceeds the shock layer equilib¬ 
rium temperature. Since the dominant ionization mechanism behind the shock front is 
electron-atom collisions (Ref. 1) that are governed by the free electron temperature, this 
enhancement of T e accelerates ionization, shortens the chemical nonequilibrium zone to 
about 0.4 cm, and decreases the overall shock layer thickness. It should be noted that 
this result shows the strong sensitivity of the overall solution to the electron temperature 
model at such trajectory points. 

Radiative heating results obtained with the present model are shown on Table I. Since 
at these conditions the radiative transfer should be dominated by atomic processes and 
since for the predicted shock layer temperatures nitrogen should be a reasonable model 
for air (Ref. 9), these results, which have been obtained utilizing only atomic processes 
corrected for LTNE effects, should be appropriate. As can be seen, the predicted stagnation 
point radiative heat transfer for the case without any radiation gasdynamic coupling is 1691 
watts/cm 2 , which is in remarkable agreement with the corresponding prediction of Ref. 6. 
It should be noted, however, that the present results indicate that most of the radiation 
originates from the high temperature equilibrium portion of the shock zone and not from 
the nonequilibrium part as postulated in Ref. 6. In the post-shock region, chemical 
nonequilibrium induces local thermodynamic nonequilibrium and depopulates the excited 
states rapidly via ionization, with the result that very little radiation originates in the 
nonequilibrium region. 

Moreover, the radiation coupled results for this case indicate significant radiation 
cooling is present. This cooling is evidenced not only by the decrease in radiative heating 
to 1039 watts/cm 2 but also by the temperature and ionization profiles on Figure 4. These 
figures compare the uncoupled and coupled results; and as shown by the steady decrease 
in temperature and in particular ionization throughout the equilibrium zone, radiation 
cooling for this case is significant and needs to be included in an analysis model. 

As mentioned, the RASLE prediction for this case was 970 watts/cm 2 . However, it is 
probable that the difference between this value and the present prediction is primarily due 
to the influence of reaction chemistry and the amount of absorption in the wall thermal 
layer. Since RASLE assumes equilibrium chemistry, it should predict more molecules in 
the wall layer and hence more absorption. This possibility is borne out by the fact that the 






RASLE results (Ref. 6) indicate that the wall thermal layer absorbs about 32% of the wall 
directed radiation while the present model indicates only about 15% is absorbed. Thus, it 
appears, that while most of the shock layer is in chemical equilibrium for this case, that 
chemical nonequilibrium effects may still be important and affect the radiative heating. 
In addition, it should be noted that the equilibrium chemistry formulation (i.e. forward 
and reverse rates or forward rate combined with equilibrium constant) strongly affects the 
heating results for this case. The results shown all used the forward rate combined with 
an equilibrium coefficient formulation. 

It is mentioned above that Ref. 6 indirectly criticizes the RADICAL model, stating 
that it improperly handles line radiation by smearing the lines and that as a result it should 
give incorrect answers. Unfortunately, this assertion may be the result of a misinterpreta¬ 
tion of the RADICAL radiation model and output. For convenience, RADICAL groups the 
radiative transfer into various wavelength regions and gives appropriate average values for 
these regions; and it is these values which are frequently plotted to show the variation of 
say wall radiative heating with wavelength or electron volts. An example of such a plot for 
the present case is shown on Figure 5, and at first glance it would appear that RADICAL 
does indeed smear lines to a significant extent. However, in the actual computation of 
the radiative transfer, RADICAL actually does for most line groups perform a line by fine 
integration; and the final result is actually the consequence of such a detailed calculation. 

When the results of the detailed RADICAL computation are plotted for the radiative 
flux to the stagnation point for this case, they appear as shown on Figure 6. Here a semi-log 
abscissa has been used in order to more vividly display the underlying continuum radiation 
as well as the lines. As can be seen, there are strong lines in the infrared region below 4 
ev find in the ultra-violet between 7 and 11 ev. However, in the vacuum ultra-violet above 
11 ev man y of the lines are actually absorbing part of the continuum flux as evidenced by 
the plots dropping below the continuum level. This absorption is also evident on Figure 7 
where the stagnation heating is plotted versus wavelength. Here the line absorption of the 
continuum radiation is very evident around 0.1 microns as is the underlying continuum. 

In general, the results shown on Figures 5-7 are very similar to Figures 3 and 4 in Ref. 
6. However, careful comparison indicates that the present results have significant radiation 
above 11 ev primarily due to free-bound continuum processes, while those of Ref. 6 have 
little or no flux in this region. In Ref. 6 this difference is attributed to the usage of the 
smeared band line model in RADICAL; but, as shown on Figures 6 and 7, RADICAL does 
include the lines in detail in this region, and most of the radiation above 11 ev is due to the 
continuum, not the highly absorbing lines. This fact, combined with the absence of any 
significant radiation between 4 and 6 ev (.2 to .3 microns) in the results presented in Ref. 






6, indicates that possibly Ref. 6 treated atomic continuum radiation differently. Thus, the 
seemingly good agreement between the methods may actually be serendipity and require 
further study. 

Radiation Model Comparisons 

Considering the possible sensitivity of various flowfield cases of interest to the radia¬ 
tive heating, it has been decided to compare several available radiation models and. if 
possible, evaluate the accuracy of RADICAL. For these comparisons, the following models 
have been used - (1) RADICAL, (2) NEQAIR (Ref. 10), (3) and a modified 8-step band 
model based upon Ref. 11. Unfortunately, complete spectral comparisons have not been 
possible since NEQAIR does not automatically include absorption effects while RADICAL 
and the 8-Step model do. Thus, it has been decided to exclude the highly absorbing ultra¬ 
violet region and to limit the comparisons to the, by comparison, transparent visible and 
infrared region above 2000 A. In addition, since LTNE factors are still under development, 
the comparison cases have assumed a constant temperature and pressure slab with the 
composition determined by the equilibrium calculation in RADICAL. For RADICAL and 
NEQAIR the gas has been considered to be air, while for the 8-Step model it has been 
considered equilibrium nitrogen at the same temperature and pressure. The results of the 
calculations are presented in Table II 

As can be seen, the results of all three models Me in very good agreement with the total 
variation in each case only being about 10 percent. This agreement is not really surprising 
since previous studies (Ref. 3-4) have shown that most models agree well in the visible 
and infrared. Thus, such comparisons and similar comparisons with experimental data in 
the visible and infrared arc probably not very definitive. Nevertheless, considering that 
RADICAL agrees with other models and considering that it has been extensively compared 
to experimental data over a wide range of conditions (Ref. 2 and 12), it is believed that 
RADICAL is an excellent and adequate radiation model for the present research. In other 
words, the problems with predicting AOTV flowfields are not associated with the primary 
radiation model. Instead the difficulties are a result of the sensitivity to chemistry, electron 
temperature modeling, LTNE factors, etc 

Velocity Effects at 80 km 

Results have also been obtained for three different velocities, 12, 14, and 16 km/sec, 
at an altitude of 80 km. These velocities are, depending upon the trajectory chosen, within 
the possible range of entry velocities associated with a Martian return vehicle. In all cases, 
the results are for the stagnation line of a 2.3 meter nose radius vehicle, the freestream is 
nitrogen, and only atomic radiation is considered. 




The temperature and composition profiles for the 12 km/sec case are shown on Figure 

8, and the radiative heating rates are listed in Table I. As can be seen on the figure by the 
continually decreasing temperature and the variation in the N“^ concentration, the entire 
shock layer at this flight condition is in chemical nonequilibrium. Immediately behind 
the shock front, which for the coupled case including LTNE factors is 10.7 cm from the 
wall, the electron temperature rises to a value several thousand degrees Kelvin above the 
expected equilibrium temperature and then gradually equilibrates with the heavy particle 
temperature. In the wall thermal layer, which comprises about 20% of the shock layer, 
deionization and recombination processes are important. 

By comparing the uncoupled uncorrected radiative heating for this case with the 
uncorrected coupled result, it is apparent that for the uncorrected radiatively coupled 
situation, which assumes that electronic states are populated according to a Boltzmann 
distribution, that there is significant radiation cooling. While not shown, comparison of 
uncorrected and corrected profiles indicates that this cooling occurs in the outer portion of 
the shock layer where the electron temperature is high. On the other hand, the corrected 
results, which include LTNE factors, shows only slight radiation coupling or cooling. 

Results were also obtained earlier using a version of the code which utilized both for¬ 
ward and reverse rates instead of a forward rate and an equihbrium expression. In those 
cases, the level of ionization was about 50% higher and the temperature profiles, particu¬ 
larly in the region near the shock front were different in that the electron temperature was 
higher. As a consequence, the radiative heating rates comparable to those in Table I were 
a factor of two to three higher. This difference was primarily due to the higher electron 
temperature, and again demonstrates the sensitivity of radiative heating to composition 
and electron temperature. 

Interestingly, the radiative heating value of 9.44 watts/cm 2 predicted for the situation 
including radiation coupling and LTNE effects is remarkably close to the 10.5 watts/cm 2 
previously predicted for this case (Ref. 1). The latter was obtained using an earlier version 
of the model before the shock jump and chemistry improvements were incorporated. In 
addition, it used the 8-Step radiation model instead of RADICAL. However, the old shock 
jump conditions yielded a lower heavy particle temperature behind the shock front, which 
when combined with the older chemistry model predicted a very similar electron temper¬ 
ature profile. As a consequence, the heating rates were similar. Again, the sensitivity of 
radiative heating to electron temperature is evident. 

The temperature and composition profiles for the 14 km/sec case are shown on Figure 

9, and the radiative heating rates are also listed in Table I. Since the freestream velocity 
is higher, the nonequilibrium zone behind the shock front is shorter than at 12 km/sec, 






occupying only the outer 50% of the 8.7 cm thick shock layer. Nevertheless, the flow is 
dominated by a nonequilibrium chemistry zone, composed primarily of N, N+, and e , and 
the wall thermal layer. While N 2 + peaks behind the shock front, it is very small over most 
of the shock layer. In addition, as can be seen in Table I, there is m the uncorrected case 
extensive radiative cooling. As shown on Fig. 9a, the profiles for various cases indicates 
that for the uncorrected case that the cooling occurs primarily in the nonequilibrium region 
behind the shock front where the electron temperature is highest. 

However, for the corrected cases the inclusion of LTNE nonequilibrium effects signifi¬ 
cantly decreases the radiation from the nonequilibrium portion of the shock layer since in 
those cases ionization processes deplete rapidly the excited atomic electronic states. As a 
result, while there is radiative cooling, what there is occurs occurs in the equilibrium por¬ 
tion of the shock layer between Y/Yshock of 0.2 and 0.6. In addition, the present results 
indicate that LTNE phenomena reduce the radiative heating by about 80% for this flight 
condition. 

Finally, with respect the 14 km/sec case, it should be noted that the same trend exists 
as in the 12 km/sec case with respect to the sensitivity to the old and new reaction rate 
treatments. As before, the electron temperature is higher in the nonequilibrium zone with 
the old formulation; and the radiative heating is approximately a factor of two higher. 

The temperature and composition profiles for the 16 km/sec case are shown on Figure 
10, and the radiative heating rates axe again listed in Table I. Here, the electron tempera¬ 
ture immediately behind the shock front is very high, having a peak value in the corrected 
coupled case of sightly over 20,000 K, as is the amount of ionization, which is about 60%. 
Likewise, due to the increase in velocity, the nonequilibrium zone is somewhat shorter. In 
this case it is about 30-40% of the 7.32 cm shock layer. 

As would be expected, the radiative heat transfer results have the same trend as those 
for 14 km/sec. The uncorrected cases exhibit significant radiation coupling and cooling, 
with almost all of the cooling occurring from the nonequilibrium portion of the shock layer 
through the shock front due to the very high electron temperature in that zone. This 
type of behavior is consistent with that shown in Ref. 3 in that reference s study of the 
sensitivity of the flowfield to radiation parameters. On the other hand, when LTNE effects 
are included (corrected results), the radiative cooling is significantly less. As before, it 
occurs primarily in the equilibrium portion of the shock zone. 

A graphical s umm ary of the 80 km radiative heating results is presented as Figure 10a. 
From this figure it is easy to see the tremendous effects the inclusion of LTNE corrections 
has on the total heat transfer for all three flight velocities. Also evident is the much lower 
amount of radiative coupling present in the LTNE corrected flows as compared to the 





uncorrected, flows. 

Finally, for all three flight velocities, the predicted radiative heating is significant 
compared to the convective heating; and, in the 16 km/sec case, the radiative heating 
exceeds the convective prediction by 70%. While the latter is only approximate in that 
only 51 points have been used in these calculations and that the wall has been assumed 
fully noncatalytic, it is probably reasonably accurate. Also, since it is anticipated that 
advanced heat shield materials can withstand only up to 70 watts/cm 2 without ablating, 
these results indicate that at 80 km ablative heat shields would be required on 2.3 meter 
nose radius vehicles at velocities of 14 km/sec and above. 

U = 16 km/sec Altitude = 75 km and 72 km 

In order to investigate altitude effects and to determine the difficulty of using the 
model under a situation on a “large” vehicle where much of the shock layer is in equilibrium, 
results have been obtained for the 2.3 meter radius body at 16 km/sec at 75 and 72 km. 
The composition and thermal profiles for 75 km are displayed on Figure 11, and in this 
case it can be seen that due to the lower altitude the chemical nonequilibrium zone is 
shorter occupying only about 15-20% of the shock layer. Likewise the entire shock layer is 
as a result of the higher pressure and density thinner than at 80 km, having a thickness of 
6.8 cm in the coupled corrected case. 

Further, as shown on Table I, even in the corrected case including LTNE effects there is 
significant radiative cooling. While difficult to detect on Figure 11, this cooling does affect 
the composition and temperature profiles in that the temperature is steadily decreasing in 
the “eq uili brium” zone between 20 and 80%. Also, the degree of ionization actually peaks 
at Y/Yshock of 0.85 and then due to cooling decreases by 50 percent before the effects of 
the wall thermal layer are encountered around Y/Yshock = 0.2. Finally, it should be noted 
that for this case, the radiative heating dominates the problem and exceeds the convective 
rate by a factor of almost five. 

The stagnation profiles at 72 km are shown on Figure 12; and at this altitude the 
predicted shock layer thickness is 6.6 cm, only slightly smaller than the value at 75 km. 
However, due to the increased pressure, the post-shock nonequilibrium chemical relaxation 
zone is considerably shorter at about 0.75 cm.; and the radiative heat transfer is approx- 
imately a factor of two larger at 758 watts/cm 2 . Further, the temperature profile in the 
equilibrium zone is steadily decreasing; and the ionization level peaks at the end of the 
chemical relaxation zone and then decreases due to radiative cooling throughout the rest 
of the shock layer. This radiative coupling effect can be observed in the concentration 
profiles by noting the steady decrease in [N^~] and increase in [N] from the end of the 






nonequilibrium chemistry region at Y/Yshock 0.9 to the beginning of the wall thermal 
layer around 0.15. Thus, for this case radiation coupling effects are important and do 
affect the temperature and composition of the shock layer. 

The trends in radiative heating with altitude for the above 72 km and 75 km cases 
along with the previous 80km case at a freestream velocity of 16 km/sec are shown in 
Figure 12a. Again, the importance of including the effects of coupled radiative cooling in 
the calculation of radiative heat transfer for this high speed case is evident. 

Comparison with Inviscid Equilibrium Results 

Even though a significant portion of any AOTV earth entry will be at altitudes where 
viscous and che mi cal nonequihbrium effects should be important, it is believed that it 
would be instructive to compare results obtained with the present model with inviscid 
equilibrium results. Such comparisons should indicate the validity of the present model 
and the limitations of equilibrium inviscid analysis predictions. As a result a limited 
number of cases have been computed in order to compare with the inviscid equilibrium 
tabular results of Sutton. The latter were obtained using a radiatively coupled solution of 
the inviscid flow equations at the stagnation point of a hemisphere and used RADICAL as 
the radiation model. This method has been compared extensively to ground test and flight 
measurements as described in Ref. 12. In comparing the Sutton values with the present 
results, it should be recognized that in addition to the inviscid-viscous and equilibrium- 
nonequilibrium differences, the Sutton results were obtained for air while the present values 
assumed a freestream of nitrogen and only include atomic radiation at this point. 

In this comparison effort, four cases have arbitrarily been selected; and these are 
listed on Table III along with the inviscid equilibrium results. The corresponding viscous 
nonequilibrium heat transfer results axe given in Table I, and the shock layer profiles are 
presented in Figures 12 - 15. 

In general, the heating predictions from the two methods are of the same order of 
magnitude; and the shock standoff distances are similar. However, there are interesting 
differences. First, the shock standoff distances from the VSL nonequilibrium solutions 
are usually less than those obtained in the inviscid equilibrium cases. In the viscous 
nonequilibrium situation, the nonequilibrium zone behind the shock front has a lower than 
equilibrium density, which would tend to cause the shock layer to be thicker than in the 
equilibrium case. On the other hand, the wall thermal layer has a very high density due to 
the assumed cool wall temperature, which would case the shock layer to be thinner. Ap¬ 
parently for the cases considered, the effects are counterbalancing with the result that the 
viscous nonequilibrium shock layer thickness is slightly less than the inviscid equilibrium 






result. 

The second interesting point is that in all cases the nonequilibrium viscous radiative 
heat transfer is less than the corresponding equilibrium inviscid value. In the nonequilib- 
rium case, LTNE phenomena significantly reduces the radiation originating m the chemical 
nonequilibrium region behind the shock front and the cool temperatures in the wall ther¬ 
mal layer reduce radiation from that zone. Since the shock layer thicknesses are similar, 
these effects reduce the radiation heat transfer prediction to values below the equilibrium 
results. 

In spite of the differences between the viscous nonequilibrium and inviscid predictions, 
it appears that the present results are reasonable and demonstrate the importance of 
including viscous, and chemical and radiative nonequilibrium effects in the AOTV flight 
regime. 


EXPECTED RESULTS 

In addition to the above results, the following additions and/or modifications to the 
flowfield code and radiation code currently being considered and are expected to be in¬ 
cluded in the final paper. First, the present approximation for the electron temperature 
will be replaced by a complete electron-electronic energy equation fully consistent with the 
approximations inherent to the VSL method. Since the added terms will be of lower order 
magnitude than the current terms being included, the current solution scheme will not be 
modified, but the new terms will be added explicitly as corrections or perturbations to the 
approximate method. The authors feel that the quasi-equilibrium equation currently being 
used is an accurate approximation and do not expect the new results to differ significantly 
from the current results. 

In addition to the existing atomic LTNE correction factors currently used m the radi¬ 
ation calculations, the authors plan to have a similar set of molecular radiation correction 
factors incorporated. Past research into molecular radiation has included such factors 
which were calculated in a manner directly analogous to the atomic correction factors. 
The authors now question the accuracy of corrections factors calculated in such a manner 
given the relatively even spacing of the electronic energy levels in such molecules as N 2 , 
the current belief of many researcher in a strong coupling between free electrons and vibra¬ 
tional energy, and the effectiveness of pre-dissociation and reverse pre-dissociation in the 
depopulation and population of the excited electronic states, respectively, as compared to 
collisional dissociation. As an alternative to the previously used method, the authors are 
presently including a quasi-steady state model into the radiation calculations to calculate 
the actual individual electronic excited state populations. 





Another on going effort is being made to include a vibrational energy model, T v ‘^, into 
the VSL flowfield code. At the higher flight velocities associated with lunar and martian 
returns the inclusion of a separate vibrational equation has not been deemed necessary 
due to the fast dissociation of the diatomic species near the shock front. For the lower 
speeds associated with some orbital transfer operations and the AFE flight experiment, 
a separate vibration equation will allow more accurate prediction of the nonequilibrium, 
relaxation phenomena behind the shock. 

Also pl ann ed for inclusion in the final paper are cases for an air freestream gas mixture 
rather than just nitrogen, solutions at various FIRE 2 trajectory points, comparisons to 
the FIRE flight data, and at least one full face solutions for a AFE type configuration. 

CONCLUSIONS 

While the present viscous nonequilibrium model still needs development and improve¬ 
ment in many areas, it, even in its present form, offers several advantages over other existing 
techniques. First, it includes viscous and chemical nonequilibrium effects. Second, it is 
reasonably computationally efficient with respect to both time and resource requirements. 
Third, it utilizes a detailed radiation model, RADICAL, which accounts for the molecular 
continuum, atomic lines, and atomic continuum phenomena. Fourth, this model has been 
modified so that the effects of local thermodynamic nonequilibrium on the radiative trans¬ 
fer axe included in the computation of the atomic radiation phenomena. Finally, fifth, the 
model includes multi-temperature effects in both the nonequilibrium chemistry and radi¬ 
ation models by computing via a free electron equation model an electron temperature. 
It is believed that this approach to the computation of radiation and chemistry effects is 
applicable for those AOTV entries for which diatomic species are insignificant over most 
of the shock layer. 

Flowfield solutions obtained with this model show a number of important dependancies 
upon the approximations and formulations. The current use in this model of a single 
reaction rate and the theoretical equilibrium constant to describe each chemical reaction 
was chosen to insure the proper species concentrations at equilibrium, but had a secondary 
effect on the magnitude of radiation due to a significant change in electron densities from 
the earlier kf-k^ model. The inclusion of an electron-electronic energy equation shows 
the presence of an electron temperature overshoot in the nonequilibrium region near the 
shock. This overshoot will inturn result in a large pulse of radiant energy from that region 
unless the effects of nonequilibrium thermodynamics are included, in which case the bulk 
of the radiation is emmitted in the near equilibrium regions of the flow profiles and the 


total radiative heat transfer to the wall is much less. 





Consequently, the present study demonstrates that the prediction of high altitude 
aerocapture vehicle flowfields is strongly dependent upon the details of nonequilibnum 
chemistry, nonequilibrium radiation, and electron temperature profiles. In addition, the 
present results indicate that these phenomena are in many cases highly coupled and inter¬ 
dependent. Finally, the present results show that for many cases of interest in aerocapture 
that radiative gasdynamic coupling is significant and that this coupling is strongly in¬ 
fluenced by radiative and chemical nonequilibrium. The final paper will delineate these 
regions and 'he extent of coupling in more detail. 
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Fig. 2 — Stagnation Profiles with Radiation Coupling and LTNE Effects 
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Comparison of uncoupled and Coupled Pro-files. Both with LTNE Effects 
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Fjg. 5 — Usual Presentation of Spectral Variation of Stagnation Point Heat Transfer 
U =14.5 Km/sec, H = 65 km, Rnose = 1 m 
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APPENDIX II 


DERIVATION OF THE SPECIES CONSERVATION EQUATIONS FOR 

MOMENTUM AND ENERGY 


In Chapter II the general equation of change was 
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resented as 
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However the summation terms can be expressed In vector 
notation as 





( 1 - 3 ) 
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However since the pressure tensor Is 


[ p 'J = FA A; t J 
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and. 
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Equation (1-6) can be written 
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Likewise since v^, t, and r are independent variables 
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By expressing F as X^/mi the last term becomes 
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The various dyacLics in this equation can he rewritten using 
the relationship that 
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Hence, Equation (1-13) becomes 
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Now in general, assuming that the elastic collisions 
are binary in nature, the volumetric rate of change in 
average value of some property fa can he expressed as 


n, a <!>■• = N ' ? 4i h 
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where Aj represents the change due to encounters 
he tneen particles of type 1 and particles of type J. 
According to Chapman and Cowllng3 4 the individual terms 
of equation (I-l6) can be expressed for elastic collisions 


as 
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where prime denotes the quantity after collision. The 
parameters b and e. describe the geometry of the encounter, 
while g is the magnitude of the relative velocity between 
the two particles. The respective velocity distribution 
functions are indicated by fand fj. For the momentum 
equation the elastic portion of the collision term is 
determined from Equation (1-17) using j> ± = m^ and 
(j)^ = mjV^. The resulting expression will be designated 

as P 1;J . 

Collisions involving charged particles, however, are 
characterized by small deflections and involve many 
particles at once. Thus, for electrons and ions the 
Boltzmann binary collision idea of Equation (1-17) is no. 
strictly applicable. If, on the other hand, this fact is 
ignored and Equation (1-17) evaluated using an appropri¬ 
ate collision cross section, the result is the same as 
that obtained by using a more exact treatment such as the 
Fokker-Planck equation.^ The reason is that the impor¬ 
tant deflections are small, simultaneous, and random and 
can be treated as if they were two-body sequential colli¬ 
sions. While this approach does encounter difficulties 
for like-particle interactions such as electron 
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self-equilibration, it should be adequate for the present 
analysis, which considers relaxation between species but 
assumes self-equilibration- to be Instantaneous. Now if 
the species continuity equation is multiplied by u 
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The second term Is 
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3y using the relationships 
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this term becomes 
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As in the case of the species momentum equation the third 
and fourth terms are zero because r, and t are inde¬ 

pendent variables. The last term on the left-hand side of 
Equation (1-1) becomes 
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In this manner the species energy equation becomes 
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where 2^ is used to represent the change in energy due to 
collisions. Now by using the fact that 
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and using the specie momentum equation to eliminate Nj_Xj_ • u, 
the energy equation becomes Pan/v 
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For the energy equation, Equation (1-17) yields for elastic 
collisions 


N} Aj ft -/// fjm ; ( v/ 2 - v, 2 y txJLtJ6 


(1-33) 


This equation can be divided into two terms as follows 
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where the first term represents the rate of energy gain by 
species i due to elastic encounters between species i and j 
because of thermal motion of the particles. The second 
part represents the rate at which work is done on species i 
due to elastic interactions between species i and J because 
of directed motion cf the particles. In this work the 
first term will be denoted as j and the second as 
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where the last term represents energy change due to 
collisions. 

Then "breaking up the pressure tensor as 

[/=,•? = f.T-17 

and expressing internal energy is terms of enthalpy 

e ; = /),- -f./P: 

Equation (1-32) becomes 
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